Estimate of the primordial magnetic field helicity 
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Electroweak baryogenesis proceeds via changes in the non-Abelian Chern-Simons number. It is 
argued that these changes generate a primordial magnetic field with left-handed helicity. The 
helicity density of the primordial magnetic field today is then estimated to be given by ~ 10 2 nb where 
nj, ~ 1CP 6 /cm 3 is the present cosmological baryon number density. With certain assumptions about 



the inverse cascade we find that the field strength at recombination is ~ 10 
coherence scale ~ 0.1 pc. 



G on a comoving 



The helicity of a primordial magnetic field is of crucial 
importance in determining its subsequent evolution and 
in assessing whether the observed magnetic fields can be 
generated by amplification of the seed field by a galactic 
dynamo Q. The average helicity density of a magnetic 
field, B, in a chosen volume V is defined as: 
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where B = V x A. 

The connection between the helicity of primordial mag- 
netic fields and baryon number is arrived at by consid- 
ering the process of electroweak baryogenesis occurring 
at the time of the electroweak phase transition ||. The 
genesis of baryons requires changes in the Chern-Simons 
number 
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where, Np is the number of particle families, W^ a , 
are the SU(2) and U(l) hypercharge gauge fields, g and g' 
are SU(2) and U(l) gauge couplings, and k = 1,2, 3. 
Changes in CS are achieved via the production and dis- 
sipation of non-perturbative field configurations such as 
the electroweak sphaleron ||, or linked loops of elec- 
troweak string |^|||, or other equivalent configurations 
(for a review of electroweak strings see [§). In the cos- 
mological setting, it is believed that such configurations 
would be produced in the false vacuum phase due to the 
detailed dynamics of the electroweak phase transition, 
and would then decay in the true vacuum phase. Once 
they decay, the baryon number so produced cannot be 
washed out by the subsequent production and dissipa- 
tion of more sphalerons or equivalent configurations. 



A feature of the baryon number producing intermedi- 
ary field configurations is that they carry fluxes of non- 
Abelian magnetic fields that are twisted or linked. The 
simplest configuration to analyze is that of two linked 
electroweak strings (Fig. 1). Here the loops carry Z- 
magnetic flux and the Z magnetic lines are linked with 
each other. 

Now consider the decay of the linked Z-string config- 
uration. One path to decay is by the breaking of the 
loops of string. This occurs by the creation of monopole- 
antimonopole pairs. The segments of string then col- 
lapse. If the loops are very large, the whole process 
would be embedded in the highly conductive cosmologi- 
cal plasma in which the electromagnetic magnetic fields 
are frozen-in. Then the electromagnetic magnetic (B) 
field lines remain linked and the final helicity of the B 
field is related to that of the original Z-strings. 




FIG. 1. In the first figure, a pair of linked loops 
of Z-string are shown. They decay by the nucleation of 
monopole-antimonopole pairs (black circles in the middle fig- 
ure) which get pulled apart by the strings. Magnetic field 
(B) lines (grey curves) run between the monopole and an- 
timonopole. The whole system is embedded in a perfectly 
conducting plasma in which the B field lines are frozen-in. 
Finally the strings have decayed, leaving behind a B field 
configuration with helicity. 

Other decay channels of the strings seem possible. 
For example, the loops could collapse before the strings 
break. In general we expect that the strings will collapse 
and decay into modes of B - since this is the only mass- 
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less gauge field in the theory - and that these modes will 
retain some of the original helicity of the Z field. For an 
order of magnitude estimation of the helicity, it is suffi- 
cient to assume that the decay channel described above 
is not exceedingly improbable^]. 

An electroweak sphaleron may also be interpreted as 
a segment of string terminating on a monopole and an 
antimonopole |||j]-[ll). The helicity is non- vanishing be- 
cause the Z-string has a twist. The sphaleron is unstable 
to the untwisting of the monopole with respect to the an- 
timonopole. This will transfer helicity from the twisted 
Z-string to the B field lines going from the monopole to 
the antimonopole as the sphaleron decays. 

Having established a connection between the genesis of 
baryons and the helicity of magnetic fields, it is straight- 
forward to relate one to the other. The de-linking process 
shown in Fig. f leads to a change in the baryon number 
by Np cos(26 w ), where 9 W is the weak mixing angle (see 
if the initial Z-string configuration has linkage — f . 
The B field flux of an electroweak monopole is &a = 
(Air/e) sin 2 W , where e is the electromagnetic coupling 
(see ) . The helicity of the final configuration in Fig. 
1 is —2$^. Hence, every baryon that is produced, causes 
a change in electromagnetic helicity 
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where we have used sin 2 9 W = 0.23, Np = 3, and e 2 /47r = 
a = 1/137. Note that only changes in baryon number are 
related to changes in helicity. We further assume that the 
initial helicity density is negligible and this allows us to 
estimate the final helicity as AH. (We are assuming that 
there are no sources besides electroweak baryogenesis for 
generating magnetic helicity in the very early universe.) 
The observed number density of baryons is ~ 10~ 6 /cm 3 
|p3| . Therefore the average helicity of the primordial 
magnetic field is estimated to be ~ — 10~ 4 /cm 3 []. In 
astrophysical units this is ~ — (10 _21 G) 2 -kpc. Note that 
this is not a root-mean-squared value for the magnetic 
helicity density but a mean local value. This is because 
our universe is observed to be made of matter and essen- 
tially no antimatter - supposedly due to the CP violation 
present in particle physics. 

Once a helical magnetic field is produced there are a 
number of circumstances under which its helicity is con- 



* Numerical simulations of electroweak strings do show fi- 
nite segments of strings that evolve by the motion of the 
monopoles at the ends J/j . Also, for Abelian fields in vacuum, 
an initially helical state asymptotically evolves to a configu- 
ration with half the initial helicity ||. 

fThis is precisely the value of the magnetic helicity assumed 



served. In the MHD approximation, the magnetic field 
in Minkowski spacetime obeys the equation 
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where v is the fluid velocity and a c is the electrical con- 
ductivity of the plasma. First, if a c is infinite, dissipation 
can be ignored and the evolution of the magnetic field 
only depends on the first term on the right-hand side of 
eq. (|j). In this case, even the local helicity - defined 
by restricting V in the integral in eq. (|l|) to volumes 
bounded by field lines O] - is conserved. In a cosmo- 
logical setting, the dissipation term can be ignored as 
compared to the effects of Hubble expansion on length 
scales above a certain critical scale called the "frozen-in" 
scale: Lf ~ y/t/4Tva c , where t is the cosmic time. Mag- 
netic fields that are coherent on scales larger than L f are 
said to be frozen-in and their helicity is conserved both 
locally and globally. 

If the fluid velocities are small, it is possible that the 
first term on the right-hand side in eq. (jij) can be ignored 
and only the dissipation term is relevant. In this case the 
field dies out exponentially fast and helicity is not con- 
served. However, if the fluid velocity is not negligible, 
the evolution of the field depends on both terms. In this 
situation, Taylor |Q has argued that the local helicity 
of the field changes due to reconnections of the field lines 
and hence is not conserved; however, the global helic- 
ity which is a sum over a lot of random local changes 
is still conserved. While there is no proof of Taylor's 
conjecture, it leads to a successful explanation of the 
"reversed field pinch". We shall assume that Taylor's 
conjecture is true. The criterion for deciding if global he- 
licity is conserved is that the magnetic Reynolds number 
TZm = 4:ira c Lv for magnetic fields on a length scale L 
and characteristic fluid velocity v, should be large. With 
cr c - T/aln(l/a) ~ T/e 2 fjjjjij and L ~ l/e 2 T, the 
length scale characteristic of gauge field configurations 
such as the sphaleron, the condition TZm >> 1 gives 
v >> e 4 /47r. Whether this condition is met depends 
on the fluid dynamics during baryogenesis. Successful 
electroweak baryogenesis requires significant departures 
from thermal equilibrium which is likely to be accompa- 
nied by large fluid velocities. Therefore we will assume 
that the fluid velocities are large enough for TZm >> 1- 
As discussed above, under these circumstances, the field 
will evolve while conserving global magnetic helicity even 
though the field is not frozen-in. 

The evolution of the field after production is a central 
problem in MHD. There are several features that are ex- 
pected on theoretical grounds [ 14| , p7|JT^ ] and on the basis 
of numerical simulations Jn],^C ] . The first feature is that 
the magnetic field tangle is expected to evolve towards a 
configuration with VxB = /3B for some constant (3. Pro- 
vided the field is characterized by a single length scale, 
this state is one of "maximal helicity" - a state in which 
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the energy is minimum subject to the constraint of fixed 
global helicity. This expectation is supported by Taylor's 
analysis of the reversed field pinch. The second feature 
is that helical magnetic fields are expected to "inverse 
cascade" - that is, energy will be transferred from small 
length scales to large length scales [ pTjp^ . [l7| p0| (though 
also see [^3|). If L(t) is the coherence scale of the field, 
the existing studies at large TZm in Minkowski spacetime 
find 



L(t) 
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where the exponent £ has been determined to be 1/2 in 
numerical studies jl9 20 1, and 2/3 in analytical studies 
under various approximations. The factor L* is the ini- 
tial coherence scale of the magnetic field and i* is the 
time scale associated with the turbulence. We will take 
t„ = L*/v* where is the typical fluid velocity. The 
third feature, as decribed earlier, is that the evolution is 
expected to conserve global helicity. 

A little care is needed in applying eq. (5) to cosmology. 
The reason is that the origin of t in eq. (5) is the time at 
the start of the simulation (t = 0) while in cosmology the 
magnetic field is produced at the electroweak epoch. Fur- 
thermore, in cosmology, the flat (non-expanding) MHD 
equations can be used provided the time used is the con- 
formal time p4| , p5[ . In a radiation dominated universe, 
the conformal time r is related to the cosmic time t by: 
t = 2tlw(t 1 / 2 — t\w). The factors of t ew have been cho- 
sen to give t = at t — t ew and a(t ew ) = 1 where a(t) is 
the cosmological scale factor. Hence, in cosmology, the 
factor (t/t*)t in eq. (||) should be replaced by (r/i*)^. 
It is more transparent and approximately equivalent to 
directly use eq. (|5|) for the first Hubble time - during 
which the expansion of the universe is irrelevant to the 
evolution of the magnetic field - and then use eq. (|5|) 
with £ replaced by £/2, together with a comoving factor, 
for further evolution in the radiation dominated epoch. 

We can now evolve the magnetic field from the elec- 
troweak phase transition to the present epoch. We focus 
on the coherence scale of the field since the field strength 
can then be estimated quite easily using the conservation 
of helicity. At the electroweak scale we have seen that the 
magnetic Reynolds number is large and hence the coher- 
ence length will grow as in eq. (@) with = L*/i>* and 
L* ~ l/e 2 T ew . We find that in one Hubble time, 
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where Tp = 10 19 GeV is the Planck energy and we 
have assumed fluid velocities: v„ ~ 1. For the slower 
estimate of the inverse cascade (£ = 1/2) this gives 
L ew ~ 1Q 8 /T ew ~ I0~ 8 t ew which is greater than the 
frozen-in scale (~ 10 7 /T ew ). For the faster inverse cas- 
cade (£ = 2/3), the result is L ew ~ 10 n /T ew . Hence 



the magnetic field becomes coherent on a scale larger 
than the frozen-in scale within one Hubble expansion af- 
ter production at the electroweak epoch. In this period 
we also expect the field to evolve towards maximal he- 
licity with the dissipation of the non-helical component 
p6[ . Once the field is maximally helical, further dissi- 
pation does not occur because such fields are force-free 

Further evolution of the magnetic field is a combi- 
nation of Hubble expansion and inverse cascade. In 
a radiation dominated universe, this gives | |l8fl L(t) = 
L{t ew )(T ew /T) 1+ £ We also check that the frozen-in scale 
Lf grows as (T ew /T) 3 / 2 . Hence the coherence scale re- 
mains greater than the frozen-in scale for both values of 
£ and the helicity of the field continues to be conserved. 

The next significant event occurs at the epoch of 
e + e~ annihilation at T ~ 0.1 MeV since the electrical 
conductivity of the plasma drops very suddenly. The 
electrical conductivity prior to the epoch is given by 
a c ~ T/e 2 (lj]]l(J while after this epoch it is given by 
cr c ~ 10 _10 m e /e 2 . This corresponds to a drop in conduc- 
tivity by a factor of 10~ 9 and an increase of the frozen-in 
scale by ~ 10 4 . The coherence scale just prior to e + e~ 
annihilation is: 
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where we have denoted quantities at e + e~ annihilation 
by the subscript ee and t ee — denotes the time just prior 
to annihilation. Inserting numbers gives L(t ee —) = 
WLf(tee-) for £ = 1/2 and L(t ee -) = 10 4 L/(i ee -) for 
£ = 2/3. Furthermore, at this stage the Reynolds num- 
ber is TZm — 4:wa c Lv = (L/Lf) 2 >> 1 where we have 
made use of the definition of Lf and have estimated the 
fluid velocity as due to cosmological expansion at the 
scale L: v ~ L/t. After e + e~ annihilation, however, 
L/L f - 10~ 3 if £ = 1/2 and L/L f ~ 1 if £= 2/3. In the 
first case, the magnetic Reynolds number is small and 
the fields are coherent on scales smaller than the frozen- 
in scale. Therefore with ^ = 1/2 we expect the field to 
get dissipated. Only the Fourier modes of the magnetic 
field larger than the frozen-in scale can survive. In the 
second case, the magnetic Reynolds number may or may 
not be significantly larger than one and the coherence 
length of the field is comparable to the frozen-in scale. 
Therefore the coherence scale of the field will grow with 
the Hubble expansion and there may or may not be fur- 
ther inverse cascade. From now on we will only consider 
the value £ = 2/3 and give estimates of the coherence 
scale of the field both with and without inverse cascade 
in the post e + e~ annihilation universe. 

The coherence scale at the recombination epoch can 
now be estimated as the scale at the electroweak epoch 
multiplied by the corresponding Hubble expansion factor 
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and the inverse cascade factor: 
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where L ew ~ 10~ 5 cms, £ = without inverse cascade 
and £ = 2/3 with inverse cascade, and the last factor 
takes into account the evolution in the matter-dominated 
era from the epoch of matter-radiation equality (T eq ~ 1 
cV) to the epoch of recombination (T rec ~ 0.1 eV). This 
gives L rec ~ 10 7 cms without inverse cascade and L rec ~ 
10 14 cms with inverse cascade. 

The strength of the field can be estimated by using the 
conservation of helicity. At recombination the magnitude 
of helicity density is given by ~ 10 2 nb rec ~ 10 5 /cm 3 or 
- (10" 13 G) 2 -(10 14 cms). (The baryon density at recom- 
bination is z 3 ec higher than that today where z rec ~ 10 3 
is the cosmological redshift at recombination.) So the 
maximally helical, primordial magnetic field with comov- 
ing coherence scale ~ 0.1 pc has a field strength ~ 10 
G at recombination. 

The present scenario also allows us to estimate the net 
helicity of the galactic magnetic field. Both magnetic 
helicity and baryon number are conserved during galaxy 
formation when the protogalactic cloud collapses. Dur- 
ing this collapse, the helicity and baryon densities in- 
crease due to the decrease in the cloud volume (V in eq. 
(Q)) but the ratio of the helicity density to the baryon 
density remains unchanged. Hence the helicity density 
of the galactic magnetic field can be estimated to be 

W 2 n b>gai 10/cm 3 which is (10- 19 G) 2 -kpc. 

Subsequent evolution of the magnetic field during galaxy 
formation, including the amplification by a turbulent dy- 
namo, does not significantly change the global helicity 
p7| and we expect this estimate to hold even today. Fur- 
thermore, we have earlier noted that the magnetic he- 
licity density has the same sign everywhere and is nega- 
tive. Hence the magnetic field in all the different galaxies 
should be left-handed^]. The left-handedness of the mag- 
netic field could also lead to a CP violating signature in 
the cosmic microwave background radiation pq| . 
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"'"This simple prediction is complicated by galactic processes 
that might generate local helicity while conserving net he- 
licity. Since then it is possible that one sign of the helicity 
may preferentially be transferred down to unobservably small 
scales. 
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